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The density of states of a 2D d-wave superconductor in the vortex state with applied 
magnetic field H in the plane is shown to exhibit fourfold oscillations as a function 
of the angle of the field with respect to the crystal axes. We find further that the 
frequency dependence of the density of states and the temperature dependence of 
transport coefficients obey different power laws, thus leading to different magnetic 
scaling functions, for field in the nodal and anti-nodal direction. We discuss the 
consequences of this anisotropy for measurements of the specific heat. 



Some of the first results supporting the d-wave symmetry of the energy 
gap in the hole-doped high-T^ cuprates came from measurements of the spe- 
cific heat in the vortex state a. This has triggered a renewed interest in the 
properties of unconventional superconductors in the mixed state nu n. 

Volovilfl showed that in a magnetic field H <C H c i the density of states 
(DOS) of a superconductor with lines of nodes is dominated by the extended 
quasiparticle states, rather than by the bound states in the vortex cores as 
for an s-wave gap. In a semiclassical treatment, the energy of a quasiparticle 
with momentum k at position r is shifted by <5wk(r) = v s • k, where v s (r) 
can be approximated by the velocity field around a single vortex. Physical 
quantities then depend on r and have to be averaged over a unit cell of the 
vortex lattice. This method has been successful in describing the thermal and 
transport properties of the vortex state with the field H perpendicular to the 
Cu0 2 layers fflil. 

Recently we have argued that the same approach remains qualitatively 
valid in YBa2Cu307-a and other not too anisotropic materials for H parallel 
to the layers and have computed the D09lJ. We assume an order parameter 
Ak = Ao cos 2(f> over a cylindrical Fermi surface parameterized by the angle <f>. 
The field H is applied in the a-b plane, at an angle a to the x-axis. The flow 
field of a vortex, v s , is elliptical due to the anisotropy of the penetration depth- 
e = A c /A a b, but after the rescaling of the c-axis v s = h(3/2mr is isotropicli3 
Here r is the distance from the center of the vortex, and /3 is a unit vector 
along the current. Approximating the unit cell of the vortex lattice by a circle 
of radius R — a -1 y/ &q/itH, where a ~ 1 is a geometric constant, we obtain 
<Wk(r) = Eh sin/3 sin(</> — ot)/p, where p = v/R and a typical magnetic energy 
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is defined as Eh = av* \JnH / '$o/2. In London theory v* = Vf/s/e, where vj 
is the Fermi velocity in the plane- 

FoLundoped YBCO, e w 5.30, and the anisotropy increases with deoxyge- 
natioiO. Taking v f = 1.2 x 10 7 cm/sa we estimate E H w 7.9ay/H KT" 1 / 2 for 
the undoped samples, while for H||c £^ = {v F /v*)E H w 18.2aV# KT -1 / 2 . 
For in-plane field fek oc sin(</> — a) and the DOS AT(cj, a) depenils on the angle 
between the field and the direction of the gap nodes, givmgL3'E2l N(0, a) /No = 
2v / 2i?ffmax(| sina|, | cosa|)/(7rA ), where N is the normal state DOS. Note 
that N(0, a) cx V~H- For w, E H <C A the density of states N(u), a) ~ 
(JVi(w,a) +iV 2 (cj,a))/2, where 

jV^.a) = f & (l + ^) , ^ ^ if z = w/Ek > 1; 

I ¥3^ [C 1 + 2a;2 ) arcsinx + 3x\/l ~ x 2 \ , if x < 1, 1 j 

for i = 1,2, £i = 2?ff| sin(7r/4— a)| and £2 = cos(7r/4— In particular 



N(w,a) 
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The frequency dependence of iV(w, a) follows different power laws for the field 
along a node or an anti-node, and consequently the specific heat coefficient 
C/T, NMR relaxation time T\T and other quantities exhibit fourfold oscilla- 
tions and a T or T 2 behavior depending on the direction of the field. 

We now compute the low temperature specific heat C(T,H) as in Rem, 
and show the result in Fig. la. Here j n = ir 2 No/3, and we have used e = 7 an.' 
E H = 0.1 A . For a = 1 and A = 200K this corresponds to H ~ 6.5T. Takin; 
7„ = 20mJ/mol K 2 , the amplitude of the oscillations m.C/"j n T for H||ab at 
T = is 0.5 mJ/mol K 2 , close to a previous estimate 2 !. This amplitude is 
reduced as T increases: at T = 0.01 A ~ 2K, it is 40% of the T = value. This 
can explain why the oscillations have not been found in the one measurement 
done for two orientations of the fielca. In an orthorhombic system the induced 
s-wave component of the gap shifts the position of the DOS minimum away 
from the 7r/4 direction, and in a heavily twinned crystal, such as used in RefEJ, 
this further suppresses the amplitude of the oscillations. 

For T <C Eh, C/T varies as T and T 2 for H\\node and H|| antinode respec- 
tively. There exists a regime Eh < r< E H where the anisotropy is washed 
out, C(H\\ab)/T cx T but C(U\\c)/T ~ const. This observation can help 
resoLsje. .some of the disagreement between the specific heat data obtained in 
RefsuoB. The results of measurements both on single crystals!! with H||c, and 
on polycrystalline samplefl are well described by C/T oc Note that due to 
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large anisotropy, the supercurrents are nearly in the a-b plane for almost all ori- 
entations of the grains with respect to HE3, so that both experiments effectively 
measure C(H||c). Since the measured specific heat is a sum of the DOS depen- 
dent and "background" contributions the analysis is rather involved. Instead, 
Revaz et affl analyzed the anisotropy SC = C(H\\c) — C(H\\antinode), inter- 
preting it as a pure vortex quantity. They found SC /T temperature dependent, 
which can be understood since it becomes T— dependent for E H ^> T > Eh- 



Wempw define C/{TE H ) = N F C (X)/A , where X = T/E H is the scaling 
variablcB. In the limit X > 1 we have F C (X) = 9C(3)A + In2/2X, similar to 



Fc(X) is shown in Fig. lb for e = 5.3. For H||c thej;rossover from C/TEh ~ 
const occurs at X c ~ 0.5, which was estimated □ to be at T j\J~H = 6.5 
KT" 1 / 2 , yielding E L H w 30%/ffKT" 1 / 2 . The crossover from small to large 
X in SC(T, H), xiccurs at X a b ~ 0.15, but the predicted I< 1 behavior was 
not found abovl T/VH w 0.55KT- 1 /! which i mp i} es e h < 4\Z# KT" 1 / 2 . 
Note that the interpolation used in Reftj gives a linear correction in X , rather 
than quadratic as in Eq. (S) , to the behavior at X <C 1 , which has led to an un- 
derestimate of the crossover scale. Notice also that the crossover in SC extends 
over a decade in X. 

The experimental values for E H and Eh are within a factor of 2 of our 
estimates, and E h /Eh is then about 3 times the predicted ratio of %fe. These 
are quantitative rather than qualitative differences, given the roughness of the 
estimates and the experimental uncertainties. Also the constant a in general 
is not the same for H||a6 and H||c. Note that for Josephson coupled planes, 
the anisptropy in the the vortex lattice constants is iJ-dependent and differs 
from £3. Finally, since the elastic modjil.i of-the vortex lattice differ for H||a6 
and H||c, the lattice contributes to JCBEaM 

Qualitative agreement between our estimates and the experimental results 
suggests that the angular oscillations of the specific heat are accessible at 
T < 2K and H ~ 10 — 20T. Such measurements would be a simple bulk probe 
of gap symmetry allowing one to map out the position of the gap nodes. 
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X < 1, H\\node, 

X <C 1, H\\antinode. 
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Figure 1: a)Normalized specific heat for H||c (diamonds), H\\antinode (circles), and H\\node 
(triangles), for e = 7 (8 ~ 0.05) and Eh = O.lAo; Inset: low-T anisotropy; b) Scaling 
function for C(H||c) (diamonds), C (H\\antinode) (full circles), C(H\\node) (full triangles), 
8C (H\\antinode) (open circles), and <5C(H||node) (open triangles). 
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